The Efficiency of Defect Production in Planar Superconductors and Liquid Crystals 
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A recent experiment that sees the spontaneous creation of magnetic flux on quenching high-T c 
films has shown that earlier null results were a consequence of the lack of saturation of the Zurek- 
Kibble causal bounds against which flux density was measured. In this letter we estimate how 
efficient the production of topological charge is in planar systems, both for the aforementioned 
experiment (when flux measures topological charge) and for an earlier experiment on planar liquid 
crystals. Agreement is good. 
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There have now been several experimental tests of the 
Zurek-Kibble (ZK) causal bounds for the production of 
defects (e.g. vortices) at continuous phase transitions. 
The ZK prediction is that the separation of defects ^dejj 
at the time of their formation, satisfies a scaling law [3, 
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where tg is the rate at which the temperature is 
quenched through its critical value T c , to is the relax- 
ation time of the long wavelength modes, and £o is the 
correlation length at zero temperature. is measured 
in units of £z, the Zurek causal length, defined through 
(|T|). ^ is the minimum size of correlated domains af- 
ter the implementation of the transition^ 0. The ef- 
ficiency coefficient / < 1 that determines the extent to 
which the causal bound is saturated will depend on the 
system, but is generally assumed to be 0(1). Typically, 
for mean-field condensed matter systems dcscribable by 
Ginzburg-Landau (GL) theory, the exponent a = 1/4. 
Experiments in annular Josephson junctions^ and non- 
linear optical systems |5| have confirmed this value, where 
it is appropriate, within errors. 

In this letter we begin by examining recent experiments 
to detect Abrikosov-like vortices at a temperature quench 
in planar high-T c superconductors by measuring the mag- 
netic flux that is spontaneously produced. The first of 
these experiments was null, with no flux observed 
above the noise in the system. However, on increasing 
the quench rate by several orders of magnitude sponta- 
neous flux was produced 0, but with low efficiency. 

The experiments of Refs. and measure net mag- 
netic flux i.e. net topological charge An, where An is 
the number of vortices n + minus the number of antivor- 
tices n_. As it stands, however, Eq^ only enables us 
to determine the total number n = n + + n_ of vor- 
tices plus antivortices. To make predictions for net flux 
requires an understanding of vortex-antivortex correla- 



tions. In the absence of explicit knowledge about corre- 
lations, both Kibble 0, and Zurek addressed this 
problem for U(l) theories, with complex order param- 
eter field 4> = pe la , by assuming that the phase angle 
a makes a random walk around the perimeter of the re- 
gion through which the flux is measured, with step length 
0(£z). This is complemented by the 'geodesic rule', that 
the shortest paths are assumed to be taken in field space 
so as to minimise gradient energies. With these assump- 
tions, the variance AN in the net flux through a surface 
with perimeter C is predicted to be 
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for a of with an overall efficiency factor /, that again 
depends on the system. The relationship (J2Jl was origi- 
nally proposed by Zurek 2] for superconductors and for 
spontaneous superflow in 4 He. 

Our first goal is to show how arises from simple 
dynamics, and to estimate the efficiency factor / for su- 
perconductors. Our second goal, postponed until later, 
is to apply the same ideas to experiments 0, that 
measure £def and AN in liquid crystals. 

Care is needed in making comparisons between the two 
types of system. In particular, for flux measurement in 
superconductors, AJV of J5J) can only be part of the story. 
We also expect to see flux due to the freezing in of the 
long wavelength modes of the electromagnetic field at the 
transition 11]. In contrast to J2J), which implies strong 
vortex-antivortex correlation, this shows strong vortex- 
vortex correlation. However, for low T c and strong type- 
II behaviour the effects due directly to the magnetic field 
are negligible, as confirmed in numerical simulations [l2l 
IT3I ]. Although larger for high-T c , its contribution to the 
variance in the net flux is still expected to be relatively 
small 0] i n the experiments of 0,0, and will be ignored. 
Further, since their contribution is positive for the case 
in point, they can only lead to an effective increase in 
efficiency. 

Whereas low-T c superconductors are well-described by 
a U(l) theory, high-T c superconductors are relatively 
poorly understood. Arguably, the effective order param- 
eter of the high-r c superconductor is also a complex field 
4>. This is seen in an idealised SO(5) model of high-T c 
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superconductors with a five-component 'rotor' order pa- 
rameter field , three of whose components are the N eel 
vector of antiferromagnetism and the remaining two the 
complex field coupling to electromagnetism. The 5*0(5) 
symmetry is explicitly broken, restored only at critical 
doping. Classically, the order parameter 'rotor' has fixed 
length, forcing vortices to have antiferromagnetic cores 
[Tq . i.e. we have a generalised X — Y model (or, in the 
parlance of QFT, a non-linear cr-model). However, fluc- 
tuations soften the constraint and, unless the system is 
very close to critical doping [171. normal Abrikosov vor- 
tices with winding number n = ±1 are produced from 
the breaking of the electromagnetic U(l) for a type-II 
system. The symmetry-breaking groundstates of the the- 
ory in this case are the S 1 of the residual 'Mexican hat' 
potential. 

EqUis tested in 7] for a = 1/4. With a/2 » 1/8 
so small, there is no way that scaling can be accu- 
rately confirmed. However, the flux that is observed 
is commensurate with provided / sa 1/16, a low 
efficiency. It was noted in Ref.0 that, if the super- 
conductor were to behave as an X-Y model, for which 
a w 0.15, rather than 0.25, we would expect an order of 
magnitude more vortices and antivortices. Since we are 
hardly seeing enough for the U(l) GL case, we do not 
consider this possibility further. If, on the other hand, 
there was no correlation between vortices and antivor- 
tices, (AiV) 2 would be proportional to area A oc C 2 , 
to give AN oc C/^def oc (tQ/to)~ <7 , in contrast to J2J). 
Equivalently, with the total number of vortices and an- 
tivortices N oc A/Q e j cx (C/^def) 2 this would imply 
AN cx iV 1 / 2 . This was originally proposed in [||, but the 
new data is totally incompatible with such a possibility. 

To show how J2J is justified and that a low efficiency is 
reasonable, we need microscopic field dynamics which, in- 
evitably, not only encode causality but also discriminate 
between systems in a way that causal bounds cannot. 
Our understanding of vortex formation is that, once the 
transition temperature has been crossed, field ordering 
is implemented by the growth in the amplitudes of the 
unstable long wavelength modes. The vortices embody 
the frustration of the order parameter field as it adopts 
different groundstatc values in adjacent domains. 

In this regard, we assume that field fluctuations can 
be taken to be Gaussian until the time t, say, when the 
back-reaction can be considered to stop field growth. For 
strongly damped systems like those considered here it can 
be taken to be effectively instantaneous. Specifically, if 
the symmetry breaking scale is \<j)\ = r), the transition is 
implemented by the time t for which (|</>| 2 )t ~ rj 2 . Be- 
cause of this we have the consequences of Gaussian fluc- 
tuations imprinted on a subsequently highly non-linear 
system. 

We identify vortices by the (line) zeroes of the field at 
this time. Writing the complex order-parameter field 4> 
as (j) — (0i + i<t>2) / V%, the total density of zeroes in the 



(12) plane (intersections of vortices plus antivortices) is 

p(x) = ^[0(x)]|e jfc a j 1 (x)9 fc 2 (x)|. (3) 

where 5 2 [0(x)] = <5[</>i(x)]<5[</!>2(x)] and €jk is the anti- 
symmetric tensor. For independently Gaussian fy, the 
ensemble average total density, to be matched to QJ, 
isJUOij 
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In (0J G(r,t) = (</>(x)</)* (0)) t is the field correlation func- 
tion at time t at separation r — |x| and the deriva- 
tives are with respect to r. It should also be noted that 
d{i) = {tdef)~ 2 from ©. 

As we are considering planar systems, for simplicity 
we shall only look at the properties of field zeroes in the 
(12) plane produced by a quench. What is relevant here, 
for high-T c experiments, is not the total density (j3J), but 
rather the topological density of field zeroes of 4>(x), which 
determines net magnetic flux. This is 



p(x) = (5 2 [0(x)]e jfe 5 j ^i(x)(9 fe 02(x). 
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in which the modulus has been dropped. 

Consider a superconductor film of area A and circum- 
ference C, enclosing n+ vortices and n_ antivortices. The 
topological charge An = n + — n_ will be zero on average, 
(p) =0. However, there will be a nonzero variance 
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at time t = t when vortices form. 

Given the density p of © and the results of |3 113 1 it 
has been shown by one of us (RR) |2£| that, if g(r,t) — 
G(r,t)/G(Q,t),then 
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This equation is the pivotal equation of this paper. The 
linear dependence on C is purely a result of the forma- 
tive Gaussian fluctuations. Although it could be argued 
that the 'geodesic rule' is reasonable, when vortices have 
classical profiles, it is not obviously the case at the early 
times that are relevant here, when vortices are appear- 
ing from a tangle of line zeroes. Nonetheless, we have 
recovered the consequences of this rule without needing 
to mention geodesies. 

On the other hand, it could be argued that the use of 
(@J to count vortices is equally unreasonable, since the 
creation of classical vortices requires the non-Gaussian 
interaction between (f>i and 02- There is an answer to 
this for global U(l) theory, and our neglect of magnetic 
field fluctuations suggests that this may provide a partial 
answer here. Essentially, at time i when most of the field 
energy is in field gradients, the energy £ associated with 
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the gradient energy of the Gaussian field in a surface 
of area A takes the form 



£ 



on using the result That is, £ ps 2nrj 2 N , where 

N = (n) is the ensemble average of n = n+ + rt_, the 
total number of vortices and antivortices respectively. A 
vortex is a tube of false vacuum whose energy/length is 
e ~ 27r?7 2 , showing that, qualitatively, the energy matches 
that of a set of classical defects whose number is the num- 
ber of field zeroes (of Gaussian fields) at the time at which 
defects are formed. It follows from the Gaussian approx- 
imation that, to a first approximation, dimensionality is 
not important if we ignore the frozen fluctuations of the 
magnetic field, as we shall. For these, dimensionality is 
important, and is a determinant in demonstrating their 
smallness |14j . 

There are further qualifications. For example, zeroes 
occur on all scales, and G(0, i) is infinite without a cutoff. 
Provisionally, we introduce coarsegraining by hand, mod- 
ifying G(r, t) by damping short wavelengths / = O(£o) hi 
D = 2 spatial dimensions as 



G(r,t) = / d 2 ke lkx P(k,t)e 



In terms of the moments of P(k, t) e 
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the density d of Q is the second moment d — 
G2 /2irGi . The cutoff is irrelevant once the 'Bragg peak' 
at wavenumber ko ~ ^,def m * ne P ower nas passed to 
long wavelength modes. This has been discussed in de- 
tail elsewhere [2(],E1 by one of us (RR). 

We note that (AA) 2 requires more than the lowest 
moments of P(k,t) needed for the total density d. A 
tedious calculation |22j shows that, approximately 
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Nonetheless, as long as kP(k,t) is strongly peaked at 
fco ~ £def then, with only one length scale, we recover 
©• Eqs. (0} and (|TT)l show that the efficiency factors / 
and / of and J3J) respectively are not related directly. 
A low efficiency in producing total flux density does not 
require a low efficiency in producing net flux. 

In order to get an idea of the size of the efficiency / for 
producing net flux we need a model for g(r,t). The van- 
ishing of ((f)di4>*) requires that g'(0,t) = 0. Empirically, 
we choose 



ff (r,t)«exp(-r 2 /2f) 



(12) 



for r sufficiently small to dominate the integral, for some 
£ = 0(£z) to be determined. At times t to such 



a Gaussian form is appropriate for the diffusion of field 
ordering (independent of dimension D) that we would 
expect from a dissipative system. Inserting this in (J7J 



d 2 x (1X7^)^2^ / d 2 x{p(x)) { (8) gives 



(AN) 2 



0.059j. (13) 



If £ = £z we find / ~ 1/4, the value obtained [23j by a 
naive discretisation of the phase angle into multiples of 
7r/2 and the 'geodesic rule'. 

For diffusive behaviour we expect £ 2 oc i, independent 
of dimension. In turn, we expect t = 0(tz), where tz is 
the Zurek causal time, the first time at which the causal 
horizon is large enough to accommodate defects. Putting 
this together, in units of £0 an d £cb we would expect 



? ~ e t z /t 



(14) 



with a coefficient of proportionality 0(1). 

In mean-field tz ~ \JtQt§ 2]. Inserting this in Ijl4(l 
and then, in turn, into (|13|) . recreates 101 with a = 1/4, 
justifying the random walk assumption of ZK in units of 
£z. This is our first general result. As we noted earlier, 
EqUis tested in for a = 1/4. 

To fix the efficiency coefficient / in J5| requires the 
coefficient of proportionality in (|14|) . Explicitly, we take 
i = btz, where b > 1 is order unity. With instabili- 
ties leading to exponential growth in field amplitudes, 
we expect b to depend only logarithmically on tQ and 
the detailed properties of the system. We note that b 
only occurs in / to the power 6 -1 / 4 , and the already 
weak logarithmic dependence on the quench rate and the 
system is totally diluted. That is, the efficiency is almost 
independent of the system. 

As a demonstration we not only assume the GL theory 
for equilibrium superconductors but assume the time- 
dependent Landau- Ginzburg (TDLG) equation for the 
transition. That is, we take 



r dt 



SF 



(15) 



for the GL quadratic energy F, where r\ a is Gaussian ther- 
mal noise, whose strength is irrelevant in defining f(r,t) 
in the linear regime. This is a reliable guide to the non- 
linear theory of low-T c superconductors^. Neglecting 
thermal fluctuations, we find 0, [2JI that, in the linear 



regime, 



5 (r,i)«exp(-r 2 t /8^), 



(16) 



for small r, provided t is sufficiently large. As a result, 
g(r, t) « exp(— r 2 /86£|) at time t, whereby the efficiency 
factor / of (0 takes the value / pa 6~ 1 / 4 /6. For the linear 
equation ( | 1 51) b can be estimated from the parameters of 
the model |2l|. Values of b between 2 and 3 are plausible, 
which gives an efficiency roughly half the canonical value 
of / = 1/4, but still larger than that seen in Q- What 
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we are saying is that, although we would not defend i|15|) 
as describing the dynamics of high-T c , lower efficiencies 
than the standard guesses based on static descriptions 
are natural. This completes our discussion of high-T c 
systems. 

In the remainder of this letter, we shall take a par- 
allel look at experiments that observe vortices at the 
U(l) symmetry-breaking interface of liquid crystals at an 
isotropic- nematic (IN) transition [jj, [l(j • For uniaxal ne- 
matic liquid crystals the bulk order parameter in the ne- 
matic phase is the director vector, with order parameter 
space RP 2 — S 2 /Z2- This allows for string defects with 
11 = ±1/2 windings. The transition is first order. How- 
ever, the anchoring of the director at the two-dimensional 
IN interface forces it to lie on a cone, whereby the order 
parameter space is a circle S , and defects have winding 
number n = ±1, corresponding to U(l) breaking. It is 
suggested in 0, that the transition could evolve by 
the same growth of unstable modes that we have invoked 
for superconductors. At least, the whole formation takes 
place almost instantaneously and the structure is com- 
pletely well focussed, implying that the transition is not 
first order. We assume that the order parameter field 
is softened by fluctuations so that defects can be identi- 
fied with its zeroes. Even if there is a crossover, rather 
than a true transition, our previous analysis goes through 
unchanged provided it is not too broad |25j , which we as- 
sume. 

In Ref. 9] the authors cool a droplet of liquid crys- 
tal through its critical temperature and measure N and 
AN through square regions at the surface, where the be- 
haviour is as above. They do not attempt to confirm (JTJ 
directly, but compare the ensemble averages AN to N. 
A fit to the form 

AN = c{N) v (17) 



gives v — 0.26 ± 0.11, and a value for c of 

c = 0.76 ±0.21. (18) 

We observe that, with the same choice of correlation 
function i|12fl . the total number N of vortices and an- 
tivortices passing through a surface, area A is, from (@J, 

N = A/2ir£ 2 , (19) 

irrespective of shape (the factor of 2-7T in the denominator 
is a direct reflection of the same factor in and does 
not assume a circular surface) . For a square surface (with 
A = (C/4) 2 ) this enables us to use (|13l) to eliminate £, 
giving 

AN w 0.77(iV) 1/4 , (20) 

independently of a. 

This value is in remarkable agreement with c of l|18[) . 
and is our second result. It is to be contrasted with the 
values of c = 0.57 and c = 0.71 that are obtained under 
the simplistic assumption of random phase distribution 
in uniform arrays of triangular or square domains, re- 
spectively [jj. 

The value v — 1/4 of 120f> . already anticipated in J3J|, 
can take no other value insofar as field ordering is driven 
by early Gaussian fluctuations before non-linearity be- 
comes important. 

We thank Emil Polturak, Ariel Maniv and Rajarshi 
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ported by the COSLAB programme of the European Sci- 
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